Global monopole surrounded by quintessence-like matter 
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We present new static spherically-symmetric solutions of Einstein equations with the quintessence- 
like matter surrounding a global monopole. These new solutions of the coupling scalar-Einstein 
equations are more complicated, which depend on the parameter of equation of state —1 < w q < — |. 
A gravitating global monopole produces a gravitational field of de Sitter kind outside the core in 
addition to a solid angular deficit. In the w q — — g case, we have proved that the solution cannot exist 
since the density of quintessence-like tends to zero if w 8 — ► — ^. As a new feature, these monopoles 
have the outer horizon depending on both Goldstone field and quintessence-like. Since current 
observations constrain —1.14 < w q < —0.93, new global monopoles have interesting astrophysical 
applications. 



1. Introduction 

The phase transition in the early universe could 
have produced different kinds of topological defects 
which have some important implications in cosmology 
Point-like defects which undergo spontaneous 
symmetry breaking arise in some theories and they 
appear as monopoles. The global monopole, which has 
divergent mass in flat space-time, is one of the most 
interesting defects. The idea that monopoles ought to 
exist has proven to be remarkably durable. Barriola 
and Vilenkin [2| first researched the characteristic of the 
global monopole in curved space-time or, equivalently, 
its gravitational effects. When one considers gravity, 
the linearly divergent mass of the global monopole has 
an effect analogous to that of a deficit solid angle plus 
a tiny mass at the origin. It has been shown that 
this small gravitational potential is actually repulsive 
Liebling Q also showed that global monopoles 
without horizon only exist for the parameter of deficit 
solid angle e < 1, while for e > \/3 no static solutions 
exist at all. The configurations for 1 < e < \/3 were 
called " supermassive " monopoles. Further, Li and 
co-workers @ have proposed a new class of cold stars 
called D-stars (defect stars). One of the most important 
features of such stars, compared to Q-stars, is that the 
theory has monopole solutions when the matter field 
is absent, which makes D-stars behave very differently 
from Q-stars. Topological defects are also investigated 
in Friedmann-Robertson- Walker space-time Q. It is 
shown that the properties of global monopoles in asymp- 
totically de Sitter/anti-de Sitter (dS/AdS) space-time 
Q and Brans-Dicke theory Q are very different from 
those of ordinary monopoles. The k-field theory, in 
which the non-canonical kinetic terms are introduced in 
the Lagrangian, has been investigated in the inflation 
scenario and the cosmic coincidence problem. Another 
interesting application of k-field is topological defects, 
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called k-defects [9j. Monopoles [lfj and vortices [11 1 
of tachyon field, which are examples of k-fields coming 
from string/M-theory, have also been studied. 

The huge attractive force between a global monopole 
M and an anti-monopole M implies that the monopole 
over-production problem does not exist, because pair 
annihilation is very efficient. Barriola and Vilenkin have 
shown that the radiative lifetime of the pair is very 
short as they lose energy by Goldstone boson radiation 
|2|. No serious attempt has been made to develop 
an analytical model of the cosmological evolution of 
a global monopole, so one is limited to the numerical 
simulations of evolution (T^ |. In the c-model approxi- 
mation, the average number of monopoles per horizon 
is Nh ~ 4. The gravitational field of global monopoles 
can lead to clustering of matter, and can later evolve 
into galaxies and clusters. The scale-invariant spectrum 
of fluctuations has been given in [l^ . Moreover, one 
can numerically obtain the CMB anisotropy (5T/T) rms 
patterns [13l ] . Comparing the theoretical value to the 
observed rms fluctuation, one can find the constraint 
of parameters in the global monopole. Therefore, the 
global monopole is considered to be relevant to structure 
formation in the early universe. 

On the other hand, current observations, such as 
CMB (Cosmic Microwave Background) anisotropy, 
SNela (Supernovae type la) and large scale structure, 
converge on the fact that a spatially homogeneous and 
gravitationally repulsive energy component, referred as 
dark energy, accounts for about 70 % of the energy 
density of universe. Some heuristic models that roughly 
describe the observable consequences of dark energy 
were proposed in recent years, a number of them 
stemming from fundamental physics [l4| and other 
being purely phenomenological [l5|. One possibility is 
that the Universe is permeated by an energy density, 
constant in time and uniform in space. In this case, 
the ratio w of pressure to energy density (the equa- 
tion of state) is also a constant at all time. Another 
possibility is that the dark energy is some kind of 
dynamical fluid, so that the ratio is a function w(a), 
where a is the scale factor of universe. Whatever, w 
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should be in the range of— 1 < w < — ^, which is the 
requirement of cosmic acceleration. Furthermore, the 
Einstein equations for the static spherically-symmetric 
quintessence-like surrounding a black hole[l6| have been 
investigated. Kiselev gave the general expression for the 
static spherically-symmetric energy-momentum tensor 
of quintessence- like matter [16(. 

In this paper, we investigated new static spherically- 
symmetric solutions of Einstein equations with the 
quintessence-like matter surrounding a global monopolc. 
These new solutions of the coupling scalar-Einstein 
equations are more complicated, which depend on the 
parameter of equation of state —1 < w q < — g. We 
show that the gravitating global monopolc produces 
a gravitational field of de Sitter kind outside the core 
in addition to a solid angular deficit. Furthermore, 
w q = — | solution cannot exist because the density of 
quintessence-like tends to zero a.s w q — ► — The new 
feature is the appearance of outer horizon for the case of 
quintessence-like matter surrounding a global monopolc. 
Especially, we studied the w q = — | case in detail. 

2. Equations of motion 

We shall work within a particular model in unit 
c = 1, where a global 0(3) symmetry is broken down to 
U(l). The Lagrangian density is 

C = \ g ^d^ a dA a - y(0> a - <to 2 ) 2 . (i) 

where <f> a is triplet of scalar fields, isovector index a — 
1,2,3. The hedgehog configuration describing a global 
monopole is 

4> a = cr /(f)^r, with x a x a = r 2 . (2) 



so that we shall actually have a monopole solution if / — > 
1 at spatial infinity and / — * near the origin. 

The static spherically-symmetric metric can be written 

as 

ds 2 = B(f)dt 2 - A{f)d? 2 - r 2 (d6 2 + sm 2 (6)dAp 2 ). (3) 

with the usual relation between the spherical coordinates 
r,9,tp and the 'Cartesian' coordinate x a . Introducing a 
dimensionless r = ergf, from (1) and (3), we obtain the 
equation of motion for / as 

7 + (|: + i<!>']''-^-^-W = °- < 4 > 

where the prime denotes the derivative with respect to r. 

When we consider the static spherically-symmetric 
quintessence-like matter surrounding a global monopole, 
the Einstein equations can be written as 

Gp, = torG(T ltv + Tp,). (5) 

where T^ v is the energy-momentum tensor for the La- 
grangian (1), and t^v is the energy- momentum tensor of 
the quintessence-like. t^ v can be characterized by a free 
parameter w q . In the static spherically-symmetric case, 
the general expression for the energy-momentum tensor 
of quintessence-like matter is given by 

Tt * = Tr r = p q , (6) 

re 6 = 7/ = -ip,(3«;, + l). (7) 

The whole situation concerning the energy-momentum 
tensor of quintessence-like matter in static spherically- 
symmetric case is discussed in great detail [l6j . The Ein- 
stein equations (5) now are ready to be written as 
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where the small dimensionless parameter e is defined by 
^v/^G^l.OSxlO- 1 ^^^). (11) 

and another dimensionless parameter p q = p q /(iQ A . 
3. The solution in simplified models 

Before solving the coupled equations for the metric 



Eqs. (8)-(10) and the scalar field Eq. (4), let us an- 
alyze an exceedingly simplified model for the monopolc 
configuration, just to manifest the main features of the 
exact solution in a naive manner. The monopolc config- 
uration is simplified by a pure false vacuum inside the 
core, and an exactly true vacuum at the exterior 



if r < 5, 

1 if r > 5. 
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The interior solution of the Einstein equations (8)-(10) can be written as follows 



2 A 



ds 2 = [1 - -^-r 2 + l^r- 3 ™'- 1 ]^ 
12 dw„ 



dr 2 
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r 2 (d6» 2 + sin 2 6 >dtp 2 ) 



(13) 
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where po is the energy density of quintessence- like matter solved as follows 
at r = 1. The Einstein equations outside the core are 
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where e characterizes a deficit solid angle. We easily see 
that the global monopole surrounded by quintessence- like 
matter generates the outer horizon of de Sitter kind at 
r = rh if — 1 < w q < — |. From the continuity of the 
metric and its first derivatives with respect to r, we have 
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Therefore, it is possible to match an interior solution to 
an exterior global monopole solution, but only with a 
negative mass M. The solution (13)-(14) is certainly not 
an exact one for the system of equations (4) and (8)-(10), 
but is an exact solution of Eqs. (8)-(10). However, we 
need to point out that the simplified version has revealed 
the principal aspects for the coupled Einstein and scalar 
field equations. First of all, making use of Eqs. (8)-(10) 
with (6) and (7), we have 



Pg = Por 



-3(w q + l) 



(16) 



where the parameter po > corresponds to the density of 
energy of the quintessence-like matter p q being positive. 
For — 1 < w q < — i , p q is divergent at origin which looks 
surprising at first glance. There is, however, no contra- 
diction. In fact, one should consider that the energy E q 
of quintessence-like matter up to a distance R away from 
the origin 
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which is indeed finite. Secondly, the negative value (15) 
for M is not in conflict with Birkoff's theorem. As is 
known to all, Birkoff's theorem affirms that the only 
static spherically-symmetric vacuum solution of Einstein 
equations is the Schwarzschild metric, and that the pa- 
rameter M in the solution is determined by the integral 



of 7* 



t\ along the source. Indeed, when R > 8, we 



have, in the simplified model (12), 
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This quantity is really positive, and the negative con- 
stant term is just the same negative effective mass of 
Eq. (15). As Barriola and Vilenkin have indicated that 
the linear term can be understood as a kind of deficit 
solid angle, quite like the gravitational effects of other 
topological defects, such as cosmic strings and domain 
walls [2]. 

4. The solutions for coupled Einstein and 
scalar field equations 

4-1. Effective mass 

Now our attention transfers from the simplified version 
to the exact solution for the system to rigorously confirm 
the features which has been born of simplified one. For 
convenience, we choose the ansatz as follows 
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where — 1 < w q < 



|. Note that for p n — 0, the exis- 
tence of solutions without horizon is restricted by e < 1 
IH]. Eqs. (8)-(10) can be formally solved by 
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and 

£M = -rfrexp[e 2 ff*rdr]. (21) 

As r — > oo, we have -B(r) = j4 _1 (r). From Eqs. (19) and 
(20), we obtain 

M A (r) = ^oJ\L + t- + _ l) 2 + p q ]r 2 dr 

_ 4 ^ r-^°r- 3 «". (22) 

which is the effective mass function of global monopole. 
A global monopole solution should have linv-joo / = 1. 
If this convergence is quick enough in flat space, Ma{t) 
will also converge fast to a finite value. The effective mass 



of the solution is determined by the asymptotic value 
lim r _> rh Ma(t) = Ma- Numerical calculations show that 
the shape of f(r) is quite insensitive to e 2 in the range 
< e < 1, therefore the above discussion is consistent. 
In order to solve the system of equations uniquely, we 
need to introduce boundary conditions as follows 

/(0) - 0,A-\r h ) = B{r h ) = 0,A(0) - B(0) = 1. (23) 



4-2. A special example : w q — — | monopole 

As a characteristic example, we are going to in- 
vestigate w q = — | global monopole. Let us first discuss 
the asymptotic behavior at the origin. Expanding the 
functions we have 



f(r « 1) = f x r + ^e 2 f lPo r 2 + 0(r 3 ),p q (r « 1) = po^ 1 + 0(r 3 ), 
A-\r « 1) = 1 - !^° r - 1( £ 2 A 2 + 6e 2 f 2 )r 2 + 0(r 3 ), 

B(r«l) = l-^r-i t 2 AV + 0(r 3 ). (24) 



where f± is a free parameter to be determined numeri- 
I 



where the effective mass Ma = lim r _» rh M A {r) is 
determined numerically. 

To go beyond the simplified model, we find out the 
gravitational field of global monopole by numerically 
solving the combined Einstein's and Goldstone field 
equations. In all numerical calculations, we can choose 
A = 1 without loosing universality. The resulting 
solutions for f(r) and p(r) are shown in Fig. 1, where 
the parameters are chosen as e = 10 -2 , w q = — | and 
po = 0,0.1,1. It is easily found that the solution for 
the Goldstone scalar field f(r) is not much different 
from the one with po = 0. In the po = case, Eqs. 
(4) and (8)-(10) are reduced to the case previously 



cally. The asymptotic behavior (r ^> 1) is given by 



(25) 



studied in great detail in @, 0, IU . We have repeated the 
calculations and confirmed the previous results. From 
Fig. 1, we find that the shape of the curve f(r) is quite 
insensitive to the value of po in the interval < po < 1 
not only asymptotically, but also close to the origin. 
Notice that the shape of /(r) is also insensitive to the 
value of e in the interval < e < 1. However, the shape 
of f(r) is sensitive if w q ^ — | and po > 1 (see Figs. 
3-4). 

Using Eq. (14) we see that the global monopole sur- 
rounded by quintessence-like matter generates the outer 
horizon of de Sitter kind at r =Th if —1 < w q < — |. Es- 
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exterior solutions are reduced to, respectively 



ds 2 -- 

and 
ds 2 = 



e 2 po 



2 A : 



r 2 ]dt 2 



dr 2 



[1 " FU 12 ' ^ [1 - e 2 p - ^r 2 } 
-r 2 (d6 2 + sin 2 6d<p 2 ). (26) 



— — jdt 2 (27) 



[l-e 2 (l + p )- 
dr 2 



[l-e 2 (l + p ) 



2Ga a M] 



r 2 (d6 2 + sin 2 9d<p 2 ). 



FIG. 1: The functions f(r) and p(r) are plotted vs the dimen- 
sionless coordinate r for different values of po in the interval 
< pa < 1. The parameters are chosen as e = 10 -2 and 
w 9 = — |. The shape of the curve f(r) is quite insensitive to 
the value of po- However, it is not universal character for the 
cases with w q — |. 




FIG. 2: The value of the radial dimensionless coordinate, 
where A _1 (r = th) = is shown as function of po for two 
different values of e and A = 1. 



pecially, the outer horizon rh = + 



(1-e 2 ) 2 



+ 



3Ap 

for w q = — I monopole in the simplified model. In the 
realistic model, we need to have the aid of numerical 
calculation. To investigate the appearance of horizon 
in the w q = — | monopole, we have fixed e = 0.0025 
and 0.01, respectively, and studied the value of zero 
of A _1 (r), r h with A _1 (r/,) = in dependence on 
po. Our results have been shown in Fig. 2. In the 
limit p — > 0, the horizon tends to infinity. Obviously, 
is a decreasing function of po. For increasing po, 
the value of r/, tends to zero as is demonstrated in Fig. 2. 



4-3. The limit of w q 



It is worthwhile to note that w„ — 



u q — — i is a spe- 
cial case. Its asymptotic behavior must be discussed 
independently. In the simplified model, the interior and 



In the limit of r — » oo, Eq. (27) implies that the met- 
ric component |<? rr | tends to a constant value [1 — e 2 (l + 
po)] -1 - Therefore, we can rescale the definition of ra- 
dius by r — > r/ y/\g rr \ so that the sphere surface gets 
the deficit solid angle since its area becomes 4irr 2 /\g rr \ 
instead of Anr 2 . Certainly, Eq. (27) is not an exact solu- 
tion for the metric around a w q = — | global monopole. 
Unfortunately, the exterior solution Eq. (27) is not con- 
sistent with the principal features of the realistic case 
w q = — |, as we shall confirm in the text below. 

Indeed, we can rigorously prove the non-existence of 
w q = — i global monopole by reductio ad absurdum. If 
there exists an exact monopole solution with w q = — |, 
we have /(0) = 0,/'(0) 7^ 0,/"(0) = and p q > 0. In 
the neighbor of origin, the coupled Einstein and scalar 
field equations can be reduced to 
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(28) 
(29) 
(30) 



From Eq. (28), we have / = fir A (°\ where f\ is an 
integral constant. Since /'(0) ^ and /"(0) = 0, we 
give A(0) = 1 and /i = .f(0). Using this result, Eq. (29) 
can be written as 



f! + p q = o- 



(31) 



which conflicts with the condition p q > 0. 

Note that we have employed the Lagrangian density 
(1) in view of the above proof, where the Goldstone field 
and quintessence-like matter are assumed to interact 
only through their gravitational influence without direct 
interaction. Although non-existence of w q = — | global 
monopole was rigorously proved, the global monopole 
may exist if we introduce an interactive term in the 
Lagrangian density (1). 



' h 4- Asymptotic behavior 



In this subsection, 



discuss the asymptotic be- 



havior of global monopolcs for — 1 < w q < 



A. At the 
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origin, the asymptotic behavior is and 
/(r«l) = / ir + 0(r- 3 ^), 
A(r<l) = 1- iJ^-r- 3 ™*- 1 + O{r ), 

p 9 (r«l) = P0r - 3 ^- 3 + O(rT). (32) 

where 

/? = < r,-r2 <9<; i' (^) The asymptotic behavior (r 3> 1) is given by 



-6w q — 2 if — | < w 9 < - 3. 
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4-5. Numerical results 

In this subsection, we still choose A 2 = 1. As 
an example of ro, 7^ — |, the resulting solution for 
f(r) is shown in Fig. 3, where the parameters are 
chosen as e = 10~ 2 ,w 9 = — i and po = 0, 10 2 , 10 3 . 
It is easily found that the solution for the Goldstonc 
scalar field f(r) is sensitive to the value of p$. It is 
inevitable that the disparity of f(r) affects the value 
of M(r). The resulting solutions of f(r) are shown in 
Fig. 4 via different w q , where other parameters are 
chosen as £ = 0.01 and p = 10 3 . Notice that the 
shape of f(r) is sensitive to the value of w q except 
asymptotical region. Obviously, the less the value of 
w q is, the more slowly the curve of f(r) tends to unit. 
The global monopole surrounded by quintessence-like 
matter generates the horizon of de Sitter kind at r = 
if — 1 < w q < — i. In Fig. 5, the value of the radial 
dimensionless coordinate, where A~ x (r = r^) = is 
shown as a function of w q for two different values of e 
and A = 1. It is interesting to find that all the outer 
horizons tend to infinite as w 9 -> -| in figure (5b). In 
Fig. 6, we plot M (r)/e 2 for different values of po, where 
M(r) = M J 4(r)/47T(To, the dimensionless parameter 
of the effective mass. Clearly, the mass function is 
negative for all r. Especially, the shape of these curves 
is considerably varied. In Fig. 7, we display M(r)/e 2 for 
different values of w q . Obviously, the mass function is 
negative for all r, the value of which increases with w q . 
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FIG. 3: The functions f(r) are plotted vs the dimensionless 
coordinate r for different values of po- The parameters are 
chosen as e = 1CP 2 , w q — — | and A = 1. The shape of the 
curve f(r) is sensitive to the value of po except asymptotical 
region. 



5. Discussion and conclusions 

Topological defects are considered to be relevant 
to structure formation in the early universe. Global 
defects are of special interest in this context since 
they have a long-range scalar field. This leads to the 
divergent mass in flat space, but renders a strong gravi- 
tational effect when the global defects are investigated in 
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FIG. 4: The functions f(r) are plotted vs the dimensionless 



coordinate r for different values of w q . 
chosen as e = 10 _2 ,po = 10 3 and A = 
curve f(r) is sensitive to the value of w q 
region. 



The parameters are 
1. The shape of the 
except asymptotical 
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-0.85, 



FIG. 6: The mass function M(r)/e (M(r) = M A {r)/Awao) 
is plotted for different values of po- The other parameters are 
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FIG. 7: The mass function M(r)/e 2 (M(r) = M A (r)/4na Q ) 
is plotted for different values of w q . The other parameters are 
e = 10~ 2 , po = 10 3 and A = 1 




-0.46 -0.41 -0.36 

(5b) 



FIG. 5: The horizon ru is shown as a function of w q for 
different values of po, where solid line, long dash line and 
short dash line correspond to po = 0.1, 1, 10 respectively. In 
these figures, is plotted in — 1 < w q < —0.83 and —0.46 < 
w q < — |, respectively. All the outer horizons tend to infinite 
as w q — » — |. 



curved space. In the global monopole case, the linearly 
divergent mass has an effective analogous to that of a 
deficit solid angle plus a tiny mass at the origin. In 
this paper, we have got the solutions of Einstein-scalar 
coupling equations for the static spherically-symmetric 
quintessence-like matter surrounding a global monopole, 
in which the scalar field and quintessence-like matter 
are assumed to interact only through their gravitational 
influence without direct interaction. 

In a fundamental sense the stability of global monopole 
arises from the nontrivial topology of the vacuum man- 
ifold Ai. The removal of the defect entails the infinite 
cost of varying all of the field hence the defect is 
classically stable. However, the global monopole was 
created in the early universe therefore it might possibly 
form thin shell like wormhole during phase transitions. 
Rahaman and colleagues [l|| have researched the thin 
shell wormhole in the context of global monopole for 
Barriola and Vilenkin solution 0. They have obtained 
the time evolution equation of the radius a of wormhole 
throat. When the initial velocity d = 0, the radius is a 
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position of static equilibrium. They have also analyzed 
the dynamical stability of the thin shell, considering 
linearized radial perturbations around static solution of 
wormhole [l8|. By applying this method, we can obtain 
a restriction on the model parameters. In our models, 
the evolution equation of throat can be written as 

a + V(a)=0 (36) 

V(a) = V BV (a) - ^a- 3 ^- 1 (37) 

where Vbv (a) is the potential for Barriola-Vilenkin case 
[l8l |. By the linearization, the stability of equilibrium 
configurations corresponds to the condition V"(ao) > 0, 
where the prime denotes derivative with respect to a. We 
have 

V"{a) = V^ v («)-^(3^ + l)(3^ + 2)e 2 po^ 3 "'^ 3 (38) 

for global monopole surrounded by quintessence- like mat- 
ter. In the — | < w q < — i case, we have V"(a) > 
7g V (a) so that the restriction of parameters could be 
loosened; in the — 1 < w q < — | case, we have V"(a) < 
Vg' v (a) therefore the restriction becomes tighter; in the 
Wq = — | case, the restriction is the same as Barriola- 
Vilenkin case. 

Our results are mainly as follows: 

(i) The global monopole solutions are possible under 
the appropriate choice of internal parameter in the 
energy-momentum tensor of quintessence-like mat- 
ter, depending on the parameter of equation of state 



-1 < w q < -|. 

(ii) In the w q = — | case, there exists the solution of 
the static spherically-symmetric quintessence-like matter 
surrounding a black hole. But the static spherically- 
symmetric quintessence-like matter surrounding a global 
monopole cannot exist in the w q — ► — | limit, because 
p q -> if w q -> -|. 

(iii) The distinguishing feature is the appearance of 
outer horizon for the case of quintessence-like matter 
surrounding a global monopole. The numerical results 
show that the radius of horizon becomes larger and larger 
as the density of quintessence-like matter decreases. 

(iv) In ro, = — | case, Eq. (25) indicates that p 
is a higher order effect in contrast with w q ^ — | 
case. Therefore, the basic features of = — | 
monopoles are similar to ordinary global monopole if we 
take < po < 1. However, the features of new global 
monopole are considerably varied in w q 7^ — | case. Since 
current observations constrain —1.14 < w q < —0.93 
new global monopoles have interesting astrophysical 
applications. 
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